Thermo— hydrodynamics 
As a Field Theory 



Jacek Jezierski 
Department of Mathematical Methods in Physics, 
University of Warsaw, ul. Hoza 74, 00-682 Warsaw, Poland 

Jerzy Kijowski 

Institute for Theoretical Physics, Polish Academy of Sciences, 
Aleja Lotnikow 32/46, 02-668 Warsaw, Poland 

Abstract 

The field theoretical description of thermo-hydrodynamics is given. 
It is based on the duality between the physical space-time and the 
"material space-time" which we construct here. The material space 
appearing in a natural way in the canonical formulation of the hydro- 
dynamics is completed with a material time playing role of the field 
potential for temperature. Both Lagrangian and Hamiltonian formu- 
lations, the canonical structure, Poisson bracket, Nother theorem and 
conservation laws are discussed. 

1 Introduction 

The simplest way to recognize the canonical structure of the electrodynam- 
ics consists in introducing the electrodynamical potentials. This way the 
homogeneous Maxwell equations become automatically satisfied and non- 
homogeneous ones can be interpreted as 2-nd order differential equations 
for the potentials. The equations can be derived from the first order, non- 
degenerate variational principle. This makes the life easy: the conservations 
laws can be deduced from the Nother theorem and the canonical structure is 
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immediately given with "Dirac delta" Poisson brackets between the poten- 
tials and their canonical momenta. The entire system becomes in a natural 
way an infinite dimensional hamiltonian system. Although it is possible to 
reduce the above structure rewriting it in terms of the electrodynamic fields 
only, without the use of potentials, but the resulting structure of the "non- 
canonical" Poisson brackets is very complicated and - for some applications 
- useless. 

The situation in hydrodynamics is very much similar. Introducing the 
material space and describing consequently the configuration of the fluid as 
a diffeomorphism between the material space and the physical space makes 
the theory formally similar to the electrodynamics. The matter continuity 
equation is automatically satisfied and the Euler equations can be interpreted 
as the second order differential equations for the field potentials. The equa- 
tions can be derived from the first order, non-degenerate variational princi- 
ple (see [6]). The corresponding hamiltonian theory is based on the canonical 
("Dirac delta") Poisson structure. 

It is very interesting that also the thermodynamics can be formulated this 
way. We introduce the potential for the temperature which can be interpreted 
as a "material time" (see [7]) and define the temperature as a ratio between 
the two different times (the material one and the physical one). This way 
the thermo-hydrodynamics becomes a lagrangian field theory which fits also 
very well into the framework of the infinite dimensional hamiltonian systems. 

Having at our disposal the two time variables, we are free to choose one of 
them as a independent variable of the theory which parameterizes the physi- 
cal events, and the other one as a field potential for the temperature. These 
two choices correspond to the "energy picture" and the "entropy picture" in 
thermodynamics. The transformation between the two pictures exchanges 
the role of the momenta canonically conjugate to the two "time variables", 
i.e. between the energy and the entropy. 

The above purely phenomenological approach has also a nice microscopic 
interpretation which we give at the end of the paper. Some applications to 
the non-conservative phenomena can be found in [5]. 
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2 Barotropic fluid: Lagrange picture 



We assume that the collection of all the points of the material can be orga- 
nized in a smooth, 3-dimensional differential manifold Y which we call the 
material space. The points of Y label the particles of the material. We in- 
troduce coordinates y a (a = 1,2,3) in the material space. Together with the 
physical time t parameterizing the time axis U they form the independent 
variables of our theory. It is useful to denote y° :— t and to consider coordi- 
nates (y a ), with a = 0, 1,2,3, in the space Y x U of independent variables 
of the theory. 

To describe different types of materials we have to equip the material 
space Y with different geometric structures. For example, the elastodynam- 
ics needs the metric (riemannian) structure in Y. For purposes of hydrody- 
namics, however, the only structure we need is the volume structure. This 
means that there is in Y a differential 3-form (a scalar density) r which mea- 
sures the quantity of the fluid (e.g. in moles) when integrated over regular 
domains of Y. In terms of coordinates y a the volume form r can be written 
as follows: 

r = h(y)dy 1 A dy 2 A dy 3 , (1) 

where the function h = h(y) is given a priori. If coordinates y a are dimen- 
sionless then h has the dimension of moles. The density h transforms simply 
when we pass from one coordinate system in Y to another: it is multiplied 
by the Jacobian of the transformation. Of course, it is always possible to 
change coordinates in such a way that the numerical value of h equals 1. 
Such coordinates will be called unimodular and we will always use them for 
the sake of simplicity. Such a formulation carries however a danger: some 
expressions have to be multiplied by h which is equal to one mole. This is 
very easy to forget. In order to have h equal to dimensionless one we will 
always choose dimension \fmole for the material space coordinates y a . This 
way the factor h can really be forgotten. 

To describe the configuration of the material at a given instant of time t 
we assign the position in the physical space X to each particle. This means 
that we have a diffeomorphism 

F : Y i — > X (2) 
describing the configuration of the material at the time t. The dynamical 
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history of the material is given by the one parameter family J 7 of configura- 
tions: 

J- : Y x U i — >X . 

If the physical space is parameterized by space coordinates x k having the 
dimension of length, then the history can be described in terms of 3 functions 
depending on four parameters: 

x k = x k {t,y a ) = x k {y a ) . 

These are field variables of the theory. The laws of hydrodynamics will be 
formulated in terms of the second order partial differential equations for the 
fields. For this purpose we need to express hydrodynamic quantities (the 
velocity and the density of the fluid) in terms of the first derivatives of the 
fields. We denote 

u dx k 

r K ■ = 

a ' dy°' 

The velocity vector v has components 

v k = 9 -^ = x k = x\. (3) 

To express the density p of the fluid we calculate the transport of the volume 
form r from the material space to the physical space. This means that in 
formula (1) we have to substitute y a as functions of x k , using the inverse 
F^ 1 : X I—)- Y of the configuration (2). Denoting 

Uk dx k ' 

where y a = y a (x k ) is the coordinate expression for the mapping we get: 

(F- l )*(hdy l A dy 2 A dy 3 ) = hdet(y a k )dx l A dx 2 A dx 3 . (4) 

We conclude that the matter density (moles per volume!) is given by the 
determinant: 

p = hdet(y a k ) = det(y a k ) , (5) 

where the matrix (y a k ) is considered not as a primary quantity but as a 
nonlinear function of the matrix (x k a ), namely its inverse. The inverse of the 
density is equal to the molar volume V of the fluid: 

V = = det {x\) (6) 



or, in terms of the volume form: 



F^dx 1 A dx 2 A dx 3 ) = Vdy 1 A dy 2 A dy 3 . (7) 

The laws of hydrodynamics will be derived from the variational formula 
with the field lagrangian equal to the difference between the kinetic and the 
potential energy: 

C = ^Mv 2 - e(V) (8) 

(if the coordinates y a are not unimodular then the above expression has to 
be multiplied by h). The constant M denotes the molar mass of the fluid (it 
could also be ^/-dependent for a nonhomogeneous fluid). The internal energy 
e = e(V) is given for the given fluid. It determines the properties of the fluid 
according to the fundamental equation: 

de(V) = -pdV , (9) 

where p denotes the pressure. Due to equations (3) and (6) the lagrangian 
density (8) can be considered as a function of field variables and their first 
derivatives together with the independent variables of the theory: 

C = C(x\x k a -y a ) . (10) 

We will show that the Euler-Lagrange equations of the above field theory 

dC 3 dC 



dx k dy a dx k a 







are equivalent to the equations of hydrodynamics. It is useful for this purpose 
to analyse the canonical structure of our theory. Introduce the momenta 
canonically conjugate to field variables: 

This way the Euler-Lagrange equations generated by the lagrangian C can 
be written as follows: 

5 ;n^|§. (12) 



dy a dx k 
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Equations (11) together with (12) can be interpreted as the generating for- 
mula for the symplectic relation (see [10] and [8]): 

Using the explicit form (8) of the lagrangian function and the fundamental 
equation (9) we obtain the following expression for the canonical momenta: 

V\ = Mv k (14) 
and therefore V° k represents the kinetic momentum density of the fluid, 

V\=V^r=Wv\- (15) 

To prove the equation (15) we used the fact that the derivative of the de- 
terminant with respect to a matrix element is equal to the corresponding 
element of the inverse matrix times the determinant itself. 

The Euler-Lagrange equations (12) can thus be written as follows: 

Observe that 

a dp dy a dp dp 

y k 



dy a dx k dy a dx k 

Moreover, we prove that the last term of (16) vanishes identically. Indeed, 
due to (7) we have: 

= F*d(-^\dx l A dx 2 A dx 3 ) = d[(y a k —)\F*(dx l A dx 2 A dx 3 )] = 

= dKVy^—ldy 1 A dy 2 A dy 3 ] = —(Vy^dy 1 A dy 2 A dy 3 . (17) 
Equation (16) takes therefore the form of the Newton equation: 

Mpv k = , (18) 

where the force on the right-hand side is given by the gradient of the pressure. 
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The dynamical equations of the theory imply the energy and momentum 
conservation laws via the so called Nother theorem. We introduce for this 
purpose the following energy-momentum tensor: 

n = V\x% - 8%C . (19) 
To prove the Nother theorem we observe that 

d dC dC dC 

q^£(x\ x k a , y a ) = — k x\ + — Xafi + , (20) 

where by x_p we denote the second derivatives: 

f) r) 2 T k 

a(3 dx? y a) dx?dx a ' 

Combining (19), (20) and the definition (11) of the canonical momenta we 
obtain: 



^_r^+— = i—v? -—)x k +(v? -—)x k = i—v p -—)x k 

dyP a dy" { dyP k dx k> a { k dx"/ af> { dyP k dx k) a ' 

(21) 

Therefore the Euler-Lagrange equations (12) imply: 

9 ;Ti + ^- = o. (22) 



dyP a dy a 

The above equality is called in the field theory Nother theorem (see [1]). We 
notice that from 3 equations of motion (12) we get 4 equations (22). This 
means that they are not independent. Indeed equations (11) and (21) imply: 

V d y P ' a)V kV ~ W 1 dx> )X « VkX »- %f>' ' dx> )X " 

dC d p 
~~dt + dy^ ° ' 

so the components of the Nother theorem corresponding to a = 1, 2, 3 imply 
its 0-component. Due to the equality (21) they are equivalent with the field 
equations. 
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Using the explicit form of the lagrangian C we obtain the following ex- 
pressions for the components of the energy-momentum tensor: 



r° = l -Mv 2 + e (23) 
and therefore T% represents the total energy density, 

Tl = pVy\v k , (24) 
T° a = Mv k x\ , (25) 

T% = S a b (e + pV - l -Mv 2 ) . (26) 

Especially interesting is the case of the lagrangians which do not depend on 
the independent parameters of the theory. In the case of hydrodynamics this 
happens e.g. when the fluid is homogeneous (both the mass and the state 
equation are the same at different points of the material). In this case the 
lagrangian depends only on x k and x k a and we obtain: 

° T*= • (27) 



The 0-th component of the equation (27) expresses the energy conservation: 

Integrating the above equality over a finite domain D and using the Stokes 
theorem we can convert the volume integral of the last term into the surface 
integral over the boundary dD of D. The integral is equal to (minus) the 
work performed by the pressure on the boundary of D. The work performed 
at the boundary is the only reason for the change of the total energy of the 
fluid contained in D. Using the identity (17) we can also rewrite the last 
term in the above equality as follows: 



dy a dx k 
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Finally, we get the energy conservation law in the standard form: 

Using (9) the "material" components of the equation (22) take the following 
form: 

- + = + v m& = ■ (29) 

The last equality is equivalent to (18) and represents the momentum conser- 
vation law of the perfect fluid. 



3 Euler picture 

In the previous section we described the configuration of a continuous medium 
as a mapping F from its material space Y to the physical space X. Now we 
invert the role of the independent parameters and the field variables of our 
theory. To describe the same configuration of the material we will use the 
inverse mapping G = F^ 1 : 

G : X i — y Y . 

The time coordinate x° = t together with coordinates (x k ) form a system of 
spacetime coordinates (x M ) (/i = 0,1,2,3) in spacetime M :— U x X. Now 
the dynamical history of the material is given by the one parameter family 
Q of configurations: 

Q : M i — > Y . 

The laws of the hydrodynamics can be expressed as second order partial 
differential equations for the field variables y a (x fl ) which are the functions 
on the spacetime M. For this purpose we express hydrodynamic quantities 
(the velocity and the density of the fluid) in terms of the first derivatives of 
the fields. We already know the expression for the density: 

P = det(y\) . (30) 
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Moreover, 



0=-x k (t,y a (t,x 1 )) = — (t,y a (t,x l )) + x k a (t,y a (t,x l ))y a (t,x l ) , 
which implies: 

v k = ^ = ~x\y a . (31) 

The above formula expresses the velocity v k in terms of the derivatives of the 
eulerian field variables y a if we consider the matrix (x k a ) not as a primary 
quantity, like in the lagrange picture, but as a nonlinear function of (y a k ), 
namely its inverse. 

It is important to notice that the above quantities satisfy automatically 
the continuity equation and therefore we do not need to postulate it inde- 
pendently. The situation is similar as in the classical electrodynamics, where 
the fields which have been derived from the electromagnetic vector-potential 
satisfy automatically the first pair of the Maxwell equations. The second pair 
becomes the dynamical 2-order equations for the potentials. We can say that 
our field variables y a play the role of the "potentials" for the hydrodynamic 
quantities (v k ,p). The analogy with the electrodynamics goes further: there 
are gauge transformations in both theories which do not change the physical 
quantities. In the electrodynamics gauge transformations consist in adding 
a gradient to the electromagnetic vector-potential. Here, gauge transforma- 
tions are given by unimodular reparameterizations of the material space Y. 
Such a reparameterization consists only in "changing names to the particles 
of the material" and does not change the values of v k and p. 

To prove the continuity equation we first introduce the matter current j 
as a vector density defined by the formula: 

3* = e^ x y\y\y\ ■ (32) 

We have 

f = ^ vX y\y\y\ = e kmn y\y 2 mV \ = det(y\) = P . 

Moreover, we have: 

Ak a _ kiiv\ 1 2 3 a _ O/ui/A 1 2 3 a , k^uX 1 2 3 a 

j y k — e y LiiJ vU \y k — e y ^y »y \y o + e v ^y vV \y « • 
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The last term vanishes being the determinant of the 4x4 matrix which has 
necessarily 2 columns equal. Therefore 



which finally implies: 



j k y\ = -py a 



Introducing now the "four-velocity" defined as follows: 

u° := 1 , 



we have 



u k := v k 



We will now show that the continuity equation for the matter current j is 
a consequence of the fact that the volume form r is closed since its external 
differential necessarily vanishes being a 4-form in a 3-dimensional space Y. 
For this purpose we prove that the transport of r to the space-time M (such 
a transport is a differential 3-form in the 4-dimensional manifold i.e. the 
vector-density in M) is equal to j. Indeed: 

g*(r) = G^dy 1 A dy 2 A dy 3 ) = y\y\y\dx» A dx v A dx x . 
Using the identity: 

dx» A dx v A dx x = e ^ A -iLj dx° A dx 1 A dx 2 A dx 3 

8x kJ 

and definition (32) we get 

g*(r) = j K -^\dx° Adx 1 Adx 2 Adx 3 . (33) 

Finally we get the continuity equation as a consequence of (33) and the fact 
that the 3-form r is closed: 

= g*(dr) = d (f^-\dx° A dx 1 A dx 2 A dx 3 ) = ^-dx° A dx 1 A dx 2 A dx 3 . 
\ ox^ J ox» 

(34) 
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Observe that replacing r by Mr we obtain the same result for the current 
Mj. This means that the following identity is valid: 

d 



dx^ 



(Mpu^) = 



even for the non-homogeneous fluid. 

Similarly as in the Lagrange picture, the field equations can be derived 
from the variational principle. Obviously, the numerical value of the action 
integral corresponding to a given configuration has to be equal in both the 
Lagrange picture and the Euler picture. This means that the new lagrangian 
L and the old lagrangian C are related via the following formula: 

Ldx 1 A dx 2 A dx 3 = G^Cdy 1 A dy 2 A dy 3 ) = pCdx 1 A dx 2 A dx 3 . 

We conclude that the numerical value of the new lagrangian equals: 

L = p(±Mv 2 - e) , (35) 

but now we have to express it in terms of the field variables y a and their 
derivatives y a using equations (30) and (31). Again, the Euler-Lagrange 
equations together with the definition of the momenta M a canonically con- 
jugate to field variables y a can be interpreted as a generating formula for the 
symplectic relation: 

5L(y a ,y^xn = ^(8^) = + (36) 

or, equivalently: 



and 



V = (38) 



The explicit form (35) of the lagrangian L implies the following expres- 
sions for the components of M a . Using the fundamental equation (9) which 
- in terms of the density p - reads 

de = ^dp , 
P 2 
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we obtain: 

0° a = -px\Mv k , 
Q k a = - px l a [Mviv k + 5f(e + P V- -Mv 2 )] . 
Therefore, the field equations can be written as follows: 

= - xk ^ + ^ pU " MVk)]= °- (39) 

Since the matrix x k a is non-degenerate, the equations can be rewritten in 
the following, equivalent form: 

° = ^ + ^™') = ^ +M ^"- < 40 > 

where we have used the continuity equation. The operator 

— •= 



dt dx» 

is called the "substantial" derivative. It is obviously equivalent to the partial 
time derivative in the Lagrange picture (i.e. with y a being constant). The 
equations (40) are called Euler equations. They are equivalent to the field 
equations (18) in the Lagrange picture. 

As in any other field theory we can also prove the Nother theorem for the 
following energy-momentum tensor: 

n = ®"aV\ - SSL • (41) 

The theorem consists in deriving the equivalence between field equations 
and the (non-) conservation laws. Its proof is completely analogous to the 
one given in the previous section. Simple calculations lead to the identity: 

8 + £ = <£<*•.- (42) 



dx^ v dx v K dx^ a dy a 
Therefore, the field equations (38) imply the (non-) conservation laws: 

+ H = ' (43) 
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Similarly as in the Lagrange picture only 3 among them are independent (this 
is due to the fact that the continuity equation has been incorporated a priori 
into the structure of the theory). Moreover, they are completely equivalent 
to field equations (39). Indeed, equation (42) implies the following identity: 

W k+ dx k)V ~W Ua d y« )VkXhV °- 
f) r)T f) r)T 

= <a^.-aF>* = &^ + a?- (44) 

Again, the symmetry of the lagrangian with respect to the space-time 
translations: 

which is physically equivalent to the absence of the external forces, implies 
the conservation laws: 

&' = °- (46) 

Simple calculations lead to the following formulae for the components of the 
energy-momentum tensor t^ v derived from the lagrangian L: 

t° = pi^Mv 2 + e) 

t k = pv k (]-Mv 2 + e + pV) 

t\ = -pMv k 

t\ = - P 5 k - pMv k vi . 

Equation = represents therefore the energy conservation law and 

equation d fl t >1 k = expresses the momentum conservation law. 

The following identities express relations between the objects used in the 
Lagrange and Euler pictures: 



©° a = -PTI 

Q k a = -pv k T° a - P x\T b a 
o 

k 



v\ = -vt° 
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V\ = Vy a l vh\-Vy a l t\ . 

We see that in the rest frame (i.e. when v k = 0) the canonical momentum 
in one picture corresponds (up to a sign) to the energy-momentum tensor in 
the other picture and vice versa. 

4 Thermodynamics 

To describe thermal properties of the fluid we need one more potential (see 
[7] and [3]). As we will see in the sequel, the new potential r can be treated 
as a "material time". We will describe the thermo-dynamical history of an 
isoentropic flow in terms of 4 field potentials depending on 4 independent 
variables. We already know that there are 2 possible choices as far as the 
"spatial" coordinates are considered: the Lagrange and the Euler picture. 
Similarly, in both cases we may choose the physical time t or the "mate- 
rial time" r as the independent variable and the remaining parameter as 
the field potential. This way we obtain 4 possible pictures of the thermo- 
hydrodynamics. In the present paper we are going to discuss only two of 
them, keeping always the physical time together with the physical space X 
and the "material time" together with the material space Y. Remaining 
two pictures have no natural relativistic counterpart and are less interesting, 
although technically it is very easy to formulate them. 

In the present section we start with the complete Euler picture, based 
on the choice of the physical space-time M as the space of independent 
parameters. To describe the space of field potentials we add new dimension 
r to the matter space Y. This way we obtain the 4-dimensional matter 
spacetime Z with r playing the role of a "material time" (having a priori 
nothing to do with the physical time t = x°). Let us introduce coordinates 
z a (a = 0, 1, 2, 3) in the material spacetime Z putting z° = r and z a = y a . 
We will describe the history of the fluid in terms of 4 potentials z a = z a (x fl ). 
The potentials give the coordinate expression for the mapping: 

Q : M i — > Z . 

The laws of the thermo-hydrodynamics will be formulated as a system of 
2-order partial differential equations for the potentials. For this purpose we 
have to express the thermo-dynamical quantities in terms of the derivatives 
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of the potentials. We already know how to do it for the hydrodynamical 
quantities j M . For the temperature T we choose the following ansatz: 



dr_. 

dx k ' 



f3u»z° u 



(46) 



with j3 being a positive phenomeno logical constant. The microscopic inter- 
pretation of the above definition will be given in the sequel but in the present 
section we adopt a purely phenomenological point of view where r is merely 
a potential for the temperature, similarly as y a were the potentials for the 
hydrodynamic quantities j M . We will prove that the choice of the Helmholtz 
free energy f(V,T) as a potential part of the Lagrangian L, i.e. 



1 

V 



-Mv 2 - f(-,T) 
2 V ; 



P 



-Mv 2 -f(p,T) 



(47) 



leads to the field equations which are equivalent with the laws of the ther- 
mo-hydrodynamics of isoentropic flows. 

Again, the Euler-Lagrange equations together with the definition of the 
momenta M Q canonically conjugate to the field variables z a can be inter- 
preted as a generating formula for the symplectic relation: 



8L{z^z^)- . hji 
or, equivalently: 

and 



9 {Q " Jza) = { i~^ a)Sza + Q " JZ ^ (48) 



dL 
~dzF„ 



d 



A 



dL 



(49) 



(50) 



dx» a dz a 

The explicit form (47) of the lagrangian L enables us to calculate the 
components of M Q in terms of the field potentials and their derivatives. 
Using the fundamental equation: 



df = —pdV — sdT = ^-dp — sdT 



(51) 



with s being the molar entropy of the fluid, we obtain: 



M O = f3psu fl , 



(52) 
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e° a = -px k a (Mv k + (3sz° k ) , (53) 
& k a = -px l a [(M Vl + ^)v k + 5 k (f - l -Mv 2 +pV)] . 

The 0-th component of field equations gives thus the entropy conservation 
law: 

Ae-.^lM^^M. (54) 

The remaining equations read: 

9 = ~\p^ a (Mv k + Psz\)u»] - JL[p X k a (f + P V - l -Mv 2 ) = 



dx» a dx^ aV ' kJ 1 dx kLr aW 1 2' 

= ~ xk ^ + ^ pU " Mvk ^ - P< Z ^(P^) = • (55) 

The last term vanishes due to the entropy conservation and the above equa- 
tions become simply the Euler equations discussed in the previous sections. 

Similarly as in any field theory we obtain the Nother theorem relating the 
translational symmetries of the system with the energy and momentum con- 
servation laws. The energy-momentum tensor of the thermo-hydrodynamic 
system is given by the standard formula: 

t\ = eV°A - . (56) 
Using the explicit expressions for canonical momenta 0^ we obtain: 

9 Vo = PiMv 2 + $sv k z\) 
Q\y\ = -p(Mv k + psz° k ) 
9 Vo = pv\\Mv 2 + fisv^ + f + P V) 

= -p[(M Vl + Psz^v' + 5 k (f + pV- l -Mv 2 )] . 

Moreover 

V°A = /WA . 

Finally, we obtain the formula for t^ x : 

t° = p{"-Mv 2 + e) 
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t k Q = pv k (]-Mv 2 + e + pV) 

t\ = - P Mv k 
t k t = - p 5 k - P Mv k vi 
with internal energy defined by 

e := / + Ts . 

As in any other field theory we can also prove the Nother theorem for the 
energy-momentum tensor t^ x . The theorem consists in deriving the equiva- 
lence between the field equations and the (non-) conservation laws. Its proof 
is completely analogous to the one given in the previous sections. Simple 
calculations lead to the identity: 

V + £ = (£e>. (57) 



dx» v dx v dx^ a dy a 
Therefore, the field equations (50) imply the (non-) conservation laws: 

V. + £ = » (58) 



dx» v dx u 

and vice versa because matrix z a v is invertible. 

We stress that now the 4 conservation laws are independent and equiv- 
alent to the 4 field equations. In particular, the energy-momentum conser- 
vation laws imply - in our formulation - the entropy conservation. This is 
due to the fact that the latter is incorporated implicite in the structure of 
the theory via the fact that the Lagrangian L does not depend on z°. 

Again, the symmetry of the lagrangian with respect to the space-time 
translations: 

d L _ Q 

dx u 

which is physically equivalent to the absence of the external forces, implies 
the energy-momentum conservation laws: 

7^ M a = • (59) 
dx» A v ; 

In the above formulation the phenomenological constant /3 may be chosen 

arbitrarily. The choice /3 — 1 is also possible. It gives us [time x temperature] 
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for the dimension of the new potential r = z°. As we shall see in the sequel, 
it is more natural to choose the dimension of /3 equal to the temperature and 
to measure r in units of time. 



5 Hamiltonian formulation of t her mo— hydro- 
dynamics 

The hamiltonian description of the thermo-hydrodynamics can be obtained 
by the following standard, field-theoretical, Legendre transformation (see 
e.g. [1] or [8]). The generating formula (48) can be rewritten as follows: 

= Q°Jz a + Q°Jz a + ^:(® k Jz a ) = 

= 6(Q°J a ) - z a SQ° a + Q°Jz a + ^(® k Jz a ) ■ (60) 
Introducing the hamiltonian density: 

H = <d\z* -L = t° 

we get another generating formula: 

-6H = Q°Jz a - z a 5Q° a + ^(^ k jn , 

where, at given time t, the function H has to be expressed in terms of the 
"hamiltonian variables" (or "canonical variables"), i.e. in terms of the fields 
z a and of the time-like component @° Q of the momenta. It is useful to choose 
the special notation for the latter: 

vr, := &° a . 

We have therefore: 

-5H = n a 5z a - z a 57r a + A-(Q k Jz a ) . (61) 
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As we already know, the numerical value of our hamiltonian is simply: 



H = p{ l -Mv 2 + e) , 

where now we consider the internal energy e as a function of p and s. It 
contains the entire information about the thermo-hydrodynamic properties 
of the fluid due to the corresponding fundamental equation: 

de = -pdV + Tds = ^dp + Tds . (62) 

P 

In order to complete the Legendre transformation we observe that the for- 
mulae (52) and (53) imply: 

• = ? p , (63) 

p k := Mv k = . (64) 

P 

Together with the formula (30), they allow us to express the hamiltonian den- 
sity H in terms of canonical variables (z a , ix a ) and their space-like derivatives. 
The formula (61) can be rewritten as follows: 

-5H(z a , z% 7i a ) = (tt q + ^® k a )6z a + Q k Jz\ - z a 5n a , (65) 
which is obviously equivalent to the dynamical equations: 

e * = oh 

V 7 I'M 



dz\ ' 

ir _^ 

a ^dx k a ~ dz<* 

and 

z = 



d7T a 

The first two equations are usually written with the help of the variational 
derivative: 

_. 6H_ = dH_ _d_dH_ 

na 5z a ' dz* dx k dz\ ' 1 ' 
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Also the last equation can be rewritten in terms of the variational derivative 
since the hamiltonian does not depend on the derivatives of 7r a : 

The theory can be also formulated in the language of infinite-dimensional 
hamiltonian systems (see e.g. [2]). The infinite-dimensional phase space of 
the system is a functional space of 8 functions (z a ,7r a ) defined on a fixed 
domain D C X in the physical space X (of course, D can also be equal X if 
we want to describe the world filled entirely with the fluid). The dynamics 
of the system is governed by the hamiltonian H D equal to the integral of the 
hamiltonian density H: 

H D := / H(z a (x), n a (x))d 3 x , 

JD 

and the generating formula (61) becomes now: 

-SH D = I njz a - z a Sn a + / ®\5z a , (68) 

JD JdD 

where dD denotes the boundary of D and by 6^ we mean the transversal 
(with respect to dD) component of the momentum 6. The above definition 
of our infinite dimensional phase space is not complete unless we specify some 
boundary conditions which enable us to annihilate the boundary integral in 
the above formula in order to obtain the infinite dimensional analog 

-5H D = [ 7r a Sz a - z a S7T a (69) 

JD 

of the finite dimensional generating formula: 

-dH(q,p) = p a dq a - q a dp a . 

The simplest way to remove the boundary integral is to impose the bound- 
ary conditions z a \gD for the configuration variables. Physically, keeping the 
values of z a \dD constant in time is equivalent to the condition v k \dD = 0. 
Moreover, controlling the value of z°\qd (not necessarily constant in time) 
means that we control the temperature on the boundary. This means that the 
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system is kept in a thermal bath. Mathematically, the above boundary con- 
ditions mean that we consider the phase space P of functions z a which fulfill 
the Dirichlet conditions on dD. The hamiltonian dynamics in P generated 
by the hamiltonian Hp describes the mixed Dirichlet-Cauchy problem for 
the field equations. Within the space P of functions fulfilling the boundary 
conditions we have 

Sz a \ dD = 

which enables us to integrate by parts during all the calculations and to 
neglect all the boundary integrals. 

Of course, controlling the Dirichlet data is not the only way to eliminate 
the boundary integrals. Using the identity 

e\6*> = *(e V) - z a t®\ 

and defining the new hamiltonian 

F D := H D - f Q\z a 

JdD 

we can rearrange the formula (68) as follows: 

-5F D = J 7rjz a - z a 5n a - J z a 5Q\ . (70) 

Now we have to complete the definition of the phase space imposing the 
boundary conditions on 0~^|aD- The hamiltonian Fd defines the dynamics 
within the space R defined this way. Such a dynamics corresponds to the 
mixed Neumann-Cauchy problem for the field equations. There are obviously 
other ways to control the boundary conditions for our field theory. Each of 
them corresponds to a different hamiltonian but only the one corresponding 
to the Dirichlet control mode is equal to the "true" energy of the system. 

In the present paper we limit ourselves only to the discussion of the 
Dirichlet control mode. In this case the generating formula (69) is equivalent 
to the dynamical equations (66) and (67). Using equations (63), (64) and 
the explicit form of the hamiltonian density H we can easily calculate the 
right hand sides of these equations. We obtain 

# = A z o v k + T 

P Pp 
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p 

and 

which is indeed equivalent to the field equations (55). 

At the end of this section let us observe that in the rest frame (i.e. when 
v k = 0) the hamiltonian generating formula reduces to the fundamental 
equation (62). Indeed, in the rest frame we have H = pe, z a = and 
therefore equation (65) takes the following form: 

d(pe) = -(n a + ^\)5z a - Q k Jz\ + z a 5n a = 

= z°6n - ®\5z\ = ^5(Pps) + (f + P V)px k a 5y\ = 

= TS(ps) + (e-Ts + ^)5p 

P 

which is obviously equivalent to the fundamental equation (62) with both 
sides multiplied by p. 



6 Poisson bracket structure of thermo-hydro- 
dynamics 

The hamiltonian dynamics of the field theory can also be equivalently for- 
mulated in terms of the Poisson bracket between the physical observables 
i.e. functionals over the phase space P. If F and G are two such functionals, 
their Poisson bracket is defined by the standard formula (see [1]): 

s F r\=[^l ^ ^ 5 J— r?n 

l ' W Jd 5z«(x) 5ir a (x) 5z«(x) 5ir a (x) ' { } 

since the z a and n a are canonical variables i.e. their Poisson bracket is equal 
to the Dirac delta distribution: 

{z a (x),n a (y)} = 8{x-y) . 
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The Poisson bracket of any physical observable F with the hamiltonian Hp 
gives the time derivative of the observable: 



F={F,H D } (72) 
and therefore also the dynamical equations can be rewritten this way: 

z a = {z a ,H D } , 
7T a = {n a , H D } . 

The reader may easily check that the above equations are equivalent to equa- 
tions (67) and (66). 

The phase space P can be factorized with respect to the following equiv- 
alence relation: we call two elements of P equivalent if they have the same 
value of the physical parameters (v k , p, s) or, equivalently, (p k , p, s) (the latter 
parameterization will be more suitable for some calculations). Denote by Q 
the corresponding quotient space i.e. the space of the equivalence classes. It 
is interesting to notice that the Poisson bracket structure can also be factor- 
ized to Q. Indeed, the formula (71) implies the following "generalized Poisson 
brackets" for the quantities defined by (63), (64) and (30): 



2 

{Pk(x),pi(y)} = - 
P 



d 

{pk(x),p(y)} = Q^ 6 ( x -y) 

{p k (x),s(y)} = iJL.8(x-y) 
{s(x),p(y)} = . 

This means that the Poisson bracket of any two functionals which are con- 
stant on equivalence classes will also be constant on equivalence classes. Fi- 
nally, we get the following formula for the Poisson bracket of the two functions 
on Q: 

rr ^ f f ( SF 5G 5G 5F \ , . . . 
{F,G}= / - ttt-tt {Pk(x),s(y)}+ 

Jd Jd \dp k (x) 5s{y) dp k (x) 5s{y) J 



24 



( 8F 5G 5G 8F \ r . , . 

( 5F 5G 5G 5F \ r . , . 

+ few " ^M^M J = 

s(x) 5 5G six) 5G d 5F 



d p(x) Spk(x) dx k 8s(x) p(x) 5pk(x) dx k 8s(x) 
5F 8 5G 5G d 5F 

Spk(x) dx k Sp(x) Spk(x) dx h Sp(x) 

4p k (x) 5F d 5G 4p k (x) 5G d 5F 



+ 



p(x) Sp k (x) dx l Spi(x) p(x) Sp k (x) dx l Spi(x) 

(there are no surface integrals left after the integration by parts because the 
elements of the configuration space satisfy the boundary conditions). 

Since the hamiltonian is also the functional on Q, we can use the above 
"generalized Poisson structure" in order to rewrite the dynamical equations 
(72) in terms of the physical observables, without any reference to the orig- 
inal phase space P (further reduction to the barotropic case can be easily 
obtained if we limit ourselves to functions which do not depend on the en- 
tropy). Some authors (see e.g. [11], [4] or [9]) postulate the above generalized 
structure a priori. Our construction shows that in fact it is the canonical 
structure of thermo-hydrodynamics and it does not depend upon someone's 
lucky guesses. 

We want to stress that "reducing" the phase space P to Q is not a sym- 
plectic reduction typical for gauge theories. The gauge transformations in 
the latter theories are defined as "canonically conjugate to Hamiltonian con- 
straints" . Therefore, the future development of the system is determined by 
initial data up to a gauge transformation only. Fixing a gauge at any instant 
of time does not determine the trajectory unless we pass to the quotient 
phase space. 

Unlike in gauge theories, there are no Hamiltonian constraints in P and 
therefore no symplectic reduction can be performed. Formally, we can call 
gauge transformations the unimodular reparameterizations of the material 
space and the additive reparameterizations of the material time on each world 
line of the fluid separately. Indeed, such a transformation does not change 
the physical meaning of the initial data. However, once fixed the initial 
gauge, the Cauchy problem for canonical variables (z a ,7T a ) can be uniquely 
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solved. This shows that the canonical structure in P is well adapted to the 
initial value problem in the thermo-hydrodynamics and the reduction to Q 
is neither natural nor necessary. 



7 Complete Lagrange Picture 

The dynamical history of the material can also be described in terms of the 
mapping: 

J= : Z i — > M 

inverse to the mapping Q used in the Euler picture. 

In this picture the laws of the thermo-hydrodynamics will be expressed 
as the second order partial differential equations for the field variables x^(y a ) 
which are the functions on the "material" space-time Z. For this purpose 
we express thermo-hydro dynamic quantities (the temperature, the velocity 
and the density of the fluid) in terms of the first derivatives of the fields. We 
have: 

dr(t(r,y a ),y a ) 
dt 

The definition of the inverse matrix element: 

o det(A) 



(A)' 1 • (73) 



° det(*y ' 
together with equation (30), implies: 

y = (x° )- 1 det«) . (74) 



Hence, 
Moreover, 
which implies: 



det«J = ^ = -L . (75) 



v* = ^ . (76) 
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The above formulae express the velocity v k , the molar volume V and the 
temperature T in terms of derivatives of the field variables x^. 

Obviously, the numerical value of the action integral corresponding to a 
given configuration is the same in both the Lagrange picture and the Euler 
picture. This means that the old lagrangian L and the new lagrangian C are 
related via the following formula: 

Cdz° A dz 1 A dz 2 A dz 3 = ^(Ldx A dx 1 A dx 2 A dx 3 ) = 

= det( x^jLdz A dz 1 A dz 2 A dz 3 = -^=dz° A dz 1 A dz 2 A dz 3 

pT 

(unlike in the purely hydrodynamic case, when both the material and the 
physical times coincide a priori, here the 4-dimensional integrals have to be 
compared!). We conclude that the numerical value of the new lagrangian 
equals: 

£ = x° ^Mv 2 -WT)) , (77) 

but now we have to express it in terms of the field variables x M and their 
derivatives x^ a using equations (73), (74) and (76). Again, the Euler-Lagrange 
equations together with the definition of the momenta V a ^ canonically con- 
jugate to field variables x^ can be interpreted as a generating formula for the 
symplectic relation: 

^ = = (^n)^ + • (78) 

Using the explicit form of the lagrangian function and the fundamental equa- 
tion 

df = -pdV - sdT 
we obtain the following expressions for the canonical momenta: 

V\ = P -^z\-{e + V V+ l -Mv 2 ) (79) 



pa 



P v 



(80) 
(81) 
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with z a ^ being a non-linear function of x^ a , namely the element of its inverse 
matrix. 

The dynamical equations of the theory imply the energy and momen- 
tum conservation laws via the Nother theorem. The corresponding energy- 
momentum tensor equals: 

r% = - 8% .£ (82) 

Using the explicit form of the lagrangian C we obtain the following ex- 
pressions for the components of the energy-momentum tensor: 

r° = -a (83) 
and therefore — T% represents the entropy density, 

T% = , (84) 
T° a = Mv k x\ - x° a (e + P V+ \Mv 2 ) , (85) 

r a b = ^ b {f + P V- l -Mv 2 ) . (86) 

Similarly as in the previous section, the zero-zero component of the 
energy-momentum tensor plays the role of the hamiltonian of the system 
with x M and IT M := V ^ being canonical variables. To perform this Legendre 
transformation we rewrite the generating formula (78) as follows: 

-5H(x», <, n M ) = (tl, + ^n)<^ + V \ 6 < ~ ^U, , (87) 
where the hamiltonian 

H := ILpd? -C = T° = -s 

has to be expressed in terms of the canonical variables. This can be done 
if we solve equations (79) and (80) with respect to "velocities" x^ . For this 
purpose we define: 

X ■= det(x k a ) 
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and 



where by x k a we denote the purely space-like part of the matrix x fl ct and by 
z a k we denote its 3-dimensional inverse. We stress that both the scalar x and 
the vector Wk are defined by the first spacial derivatives of the configuration 
variables x^. 

Due to the following identities: 

_ V o 

X — rp Z J 

A 

z 

we may rewrite definitions (79) and (80) in terms of the quantities \ and Wk- 

Tlo=px-(e + pV+^Mv 2 ) , (88) 

il k = Mv k - pxw k ■ (89) 
Moreover, we observe that the matrix 

h 

k _ X X a _ k _ k 

x a n _ a u x a 

X o 

is equal to the 3-dimensional inverse of the matrix {y\). Indeed: 

r .k „ti \ 



X Q 



Therefore 

V = det(x\ - v k x\) = X det(5f - v k w{) = X (l - w k v k ) . (90) 
From (88), (89) and (90) we get 

MV = x(M - w k U k - pxw 2 ) 

and 

I! 2 p 2 w 2 x 2 
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Solving the above two equations with respect to e and V we can express 
them in terms of "legal" hamiltonian variables, i.e. x M (together with its first 
spacial derivatives) and 11^. 

The above equations contain only two unknowns since for given constitu- 
tive equations (i.e. for a given fluid) the pressure p is a function of e and V. 
The relation between those quantities is given by the following fundamental 
equation: 

ds(e,V) = ±de + ^dV , (91) 

or, equivalently 

}__ds 
T~de 

and 

ds ,ds x i 

Finally, the hamiltonian of the system is equal to — s(e, V) with e and V 
expressed in terms of the canonical variables. 

Let us observe that in the rest frame (i.e. when v k = 0) the hamiltonian 
generating formula reduces to the fundamental equation (91). Indeed, in the 
rest frame we have x — V, n = e and x k = 0. Therefore, the equation (87) 
reduces to 

$s = ^ a /< - x°6U = p - ^V^o) + ^Se = 

= V (sQ - V5(^ + ±5e = ^5V + ±6e . 

The equation (91) has been obtained from (62) dividing both sides by 
T. Both (62) and (91) can be interpreted as symplectic relations between 
thermo-hydrodynamic quantities characterizing the material. In the first 
case the energy e = e(s, V) is a generating function of a lagrangian sub- 
manifold in the symplectic space of parameters (s,V,T,p), equipped with 
the symplectic form dT A ds + dV A dp and with (s, V) chosen as control 
parameters. In the second case the same physical relation is described by 
the generating function s = s(e, V) in the symplectic space of parameters 
(e, V, If, |;), equipped with the symplectic form d(^) Ade + d(^) A dV and 
with (e, V) chosen as control parameters. The transformation between those 
two formulations of thermodynamics, called respectively the "energy picture" 
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and the "entropy picture" , is rather obscure. It has nothing to do with the 
Legendre transformations (e.g. the Helmholtz transformation leading from 
(62) to (51) and similar transformations relating the energy to the enthalpy 
or to the Gibbs free energy). Legendre transformations consist always in 
exchanging a control parameter (e.g. the entropy) with the corresponding 
"response parameter" (e.g. the temperature). They always refer to a fixed 
lagrange submanifold within a fixed phase space, which is described with 
respect to two different "control modes" (the comprehensive description of 
these phenomena can be found in [8]). 

Unfortunately, the transformation mixing the phase space parameters 
with the generating function itself is of different nature and we do not know 
in the literature any convincing interpretation of such a transformation. Here, 
we obtain such an interpretation in a natural way: there are two possible de- 
scriptions of the dynamics depending on whether we parameterize the phys- 
ical events by the physical time or by the material time. The transformation 
between these two parameterizations exchange also the role of the momenta 
canonically conjugate with both times (the energy and the entropy). 

8 Relativistic theory 

The relativistic version of the above theory can be found in [3]. Here we give 
the short review of the results. 

In Euler picture the relativistic thermodynamics is a field theory with four 
field potentials z a = z ^) defined on the physical space-time M equipped 
with the metric structure (not necessarily flat!). The structure of the 
material space-time remains the same as in non-relativistic case. Similarly, 
we introduce the matter current 

j := g*{r) 
which is automatically conserved since 

dj = dQ*(r) = g*(dr) = . 
In terms of coordinates the components of j are given by the formula 
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which implies the continuity equation 

d 



Decomposing now 



f = • 



where the velocity vector is normalized, we obtain the definition of the 
hydrodynamic quantities (p, u^) in terms of the first derivatives of the three 
"spatial" field potentials z k (we stress that p is now the rest frame density). 
The fourth potential is used again to define the temperature as a ratio be- 
tween the two times running along the world lines of the fluid: 

T := $u»z\ . 

We choose the rest frame free energy density as a Lagrangian of the theory: 

L = L(z^) = -pf reL (p,T) . 

There is no kinetic energy term since the relativistic free energy f re i contains 
also the rest mass contribution to the energy: 

f reL :=f + M 

(we use geometric units with the speed of light being equal to 1). The 
Euler-Lagrange equations together with the definition of the momenta M a 
canonically conjugate to field variables z a can be written as a generating 
formula for the symplectic relation: 

6L(z«, z^ xn = ^(^Jz a ) = i^a)Sz a + • 

To check that the above field theory really describes the thermo-hydrodynamics 
we calculate the components of the energy momentum tensor given again by 
the formula (56). After some calculations we obtain: 

t\ = -^det(s~) ( P eu»u u + p{5^ + u»u v )) . 
The Nother theorem now reads: 

v M n = 
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which, together with the entropy equation, is really equivalent to the equa- 
tions of thermo-hydrodynamics. 

The hamiltonian formulation and canonical structure of the theory can 
be obtained from the construction very much similar to the non-relativistic 
one. 

To describe the interaction of the hydrodynamic matter with the gravi- 
tational field it is sufficient to add the gravitational lagrangian to the above 
matter lagrangian and to take into account the gravitational degrees of free- 
dom, describing the geometry of M. Such a theory of a self-gravitating fluid 
becomes especially simple in the complete Lagrange picture. The main advan- 
tage of this picture consists in the fact that the "gravitational gauge" and the 
"hydro-thermodynamic gauge" eliminate somehow each other and the theory 
can be formulated in terms of purely physical (gauge free!) quantities. This 
way we obtain the theory which is well adapted to the numerical analysis 
and can be used e.g. to computer simulations of the geometrodynamics in 
both cosmology and astrophysics (see again [3]). 

Another advantage of our approach consists in a rather straightforward 
generalization of the hydrodynamics to the elastodynamics. The relativistic 
formulation of the latter has never obtained any satisfactory formulation. 
Here, the introduction of the riemannian structure into the material space Y 
enables us to generalize the theory in such a way that it is relativistic and 
reduces in the rest frame to the non-relativistic, nonlinear elastodynamics. 
The paper containing these results will be soon published. 

9 Microscopic interpretation of the mate- 
rial time 

We suppose that the liquid is composed of molecules with mass m. If the 
temperature of the fluid equals T, the molecules move chaotically around 
the theoretical world lines of the fluid (lines tangent to the vector field u^) 
and the mean kinetic energy of this motion with respect to the rest frame 
is equal (for the low temperature) to \¥T = ^mv 2 . Due to this motion the 
proper time t for the particles is retarded with respect to the physical time 
x° calculated along u^. For velocities v much smaller than the velocity of 
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light c this retardation can be calculated from the formula: 




^ 2 j o m 1,2 \ o o o 

--dx *(1-^)X =X ~X^. 

We identify the parameter r with the proper time retardation: 



r := x 



o 



1 



f = — T°T 



where the constant (3 = has the dimension of temperature. Hence, 

^ = T 

dx° 

similarly as in formula (46). We interpret therefore the "material time" as 
the "proper time retardation" due to the chaotic motion of the particles. 
This phenomenon enables us to construct (at least theoretically) a "radium 
thermometer" . We inject a drop of the radioactive radium into the fluid. Due 
to the chaotic motion of the particles the lifetime of the radium gets length- 
ened proportionally to the temperature of the fluid. Therefore, measuring 
the lifetime we measure the temperature. 
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